
I. INTRODUCTION 

There is evidently a great deal of activity and coqubo- 
tion in deep space, unrevealed by the naked eye or the deli- 
cate twinkling of the stars. As far as we can tell, quasars 
(QSOs) and "active galactic nuclei" (A6N) have relatively 
enormous power outputs produced in very small volumes. A 
typical quasar can produce a hundred to a thousand times the 
luminosity of a normal galaxy from a region one hundred 
thousand times smaller in size. Roughly speaking, if the 
city of Boston were a galaxy in terms of its power output and 
size, then a quasar would have the power of the entire United 
States produced in a region the size of a baseball. 

At various times, it has been suggested that new kinds 
of physical laws or phenomena must be required to explain 
these objects. But more and more evidence has accumulated 
that shows a continous range of energetic activity, starting 
with normal galaxies, going through Seyfert galaxies and 
other AGN, and joining with the less luminous quasars. 
Moreover, the hypothesized model of gas falling through a 
deep gravitational well, probably caused by a massive black 
hole, together with a net angular momentum in the bottom of 
the well, seems capable of explaining, at least qualitative- 
ly, the full range of observed phenomena, from the energetics 

to the striking, long-lived "jets" of matter observed to ema- 
* 

Based on a lecture given at the Goddard Workshop on X-ray Astronomy 
(October 1981). 
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nate from the centers of these objects. In any case# there 
does seem to be good evidence for some universality in the 
mechanisms and physical conditions. 

Among the most important issues in understanding QSOs 
and AGN are (1) the nature of the power source# (2) the radi- 
ation processes# and (3) the mechanism for formation and col- 
limation of the jets. He will discuss these issues in turn# 
giving a brief# model-independent sketch of some of the im- 
portant theoretical ideas and observations (with more em- 
phasis on our own interests)# and suggesting future work. It 
is possible that no new observations within the foreseeable 
future will satisfactorily pin down the above issues. 

We will be particularly concerned with phenomena that 
produce X-rays. The clear association of strong X-rays with 
QSOs and AGN and the rapid time variability seen in X-rays 
indicate that this region of the spectrum may contain much 
information about the conditions near the central region of 
the objects. Of course# simultaneous observations of all re- 
gions of the spectrum may provide important clues to the me- 
chanisms at work. 

For some recent reviews of this subject# see e.g. Rees 
(1977#1978,1980) # Fabian and Rees (1979)# and Bradt (1980). 

II. NATURE OF THE POWER SOURCE: ACCRETION ONTO A MASSIVE 
BLACK HOLE 

A. General Considerations 


We will tentatively adopt gas accretion onto a massive 
black hole as the "standard model" for the power source. The 
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lengch scale is then set by the Schwarzschild radius of the 
black holer ^ - 2<Stt/c z , where M is the mass of the hole. We 
will use the notation Hg i Unless the hole is rotat- 
ing near maximum angular velocity, we can expect the energy 
production to peak in the region . The luminosity. 


L, can be written as 

L= e tic 2 <•> 

* • 

where £ is the efficiency parameter and W is the mass ac- 
cretion rate. The virial temperature, T v , set by the pro- 
ton rest mass energy, is 

Ty, • io ,2 k C.2-) 

The period of a circular orbit, P, also denoted by , in 

the Newtonian approximation, is 

?*t K =* 3*lo r s Kg Cr/lOr^ C3a) 

The light travel time across this region is 

At u 3 . t/c - lO^'s (ir/or<,V (3 b) 

The total mass accreted, P^acc. f in th e active lifetime, 

n.«5 I0 7 f^e ( ( -.0*^ s-’) ^ or) . 

Both the combination of inferred active lifetimes and ob- 


1.3 b) 


served luminosities, and the sizes and velocity dispersions 
obtained from optical studies, indicate an and an M in 

f> ' hS 

the range /0 — /0 . We denote the ratio of radial vel- 

ocity, \/ r , to Kepler velocity, N/^ , by 

P =. v,/v K . CSV) 

The radial infall time, -fc r . , is then 

f ~ / 1 $ . 

Note that ^ is related to the viscosity parameter OC of the 
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"06 -model" disks (Shakura and Sunyaev 1973) by 

g = oC ('*/r >“ <Sc) 

where h is the disk thickness at radius r. In "thick accre- 
tion disks" we might expect h^r, so that . The con- 

tinuity equation gives an ion number density ,N, 

H- tlj‘ (-)"&' 0-Ant) C*Ao 

where we have normalized L in terms of the Eddington luminos- 
ity 

L Ot * trS^CKl p /<^ = 10"^ i- 1 . 

We will denote the mass of the proton and electron by m p and 
m, respectively, and <7y is the Thomson cross section. From 
equation (6) we obtain the "ion scattering depth" 



* 


N 


» ' N <* = 2 (&,) &)> ^ 


The Thomson scattering depth is 


= r (n, * n.')C^ >/ z * 



where f>_ and n + are the electron and positron number densi 
ties. 


Although Z fi might be expected to approach unity in 
sources radiating at near the Eddington limit, equation (7) , 
the large optical polarization seen in some QSOs and BL Lac 
objects (Stein, O'Dell, and Str ittroatter 1976; Angel 1978), 
together with the as yet unobserved evidence of strong Corap- 
tonization, suggests V+y, 4. I . We then obtain an upper 
limit for the mass in the central, emitting region, 


(Rees 1977) 


We* * % JZr 


rOpli 


< 4* *10^ n 


A*- fJL ' 
h Lior^i 



Since << M , most of the mass has already col- 

lapsed within the emitting region, arguing against some 
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non-black-hole models for the central power source. Of 
course , the majority of QSOs and AGN do not show strong opti- 
cal polarization. 

B. Gas flow models 

The principal requirements for efficient conversion of 
gravitational energy into radiation, £ >0*1 , are that (a) 
the viscous dissipation time scale be comparable to or short- 
er than the radial infall time scale and that (b) the cooling 
time of the gas be shorter than the radial infall time scale. 
When angular momentum is present , an accretion disk forms; 
most accretion disk models proposed satisfy the above re- 
quirements. Even when the infalling gas does not have much 
net angular momentum r a disk may form out to a radius 
fA.|0 Z fs due to the "dragging of inertial frames" by a rapid- 
ly rotating black hole (Bardeen and Petterson 1974) . 
Spherical accretion models require rapid dissipation, since 
the infalling matter is not delayed by angular momentum. 

For a review of the types of ctrdisk models, see Eardley 
et. al. (1978). These models differ in whether they are 
gas or radiation pressure dominated, optically thick or thin, 
equal or unequal temperatures for electrons and ions, and in 
the source of photons; they have in common the assumption of 
a spatially constant value for the viscosity parameter 0 L and 
the assumption of a thin disk ln« * . Recently, Lynden-Bell 
(1978) and Jaroszynski, Abramowicz, and Paczynski (1980) have 
considered geometricaly thick disks, but without including 
the effects of viscosity. These thick-disk models are char- 
acterized by an unspecified free function, l(r), the 
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non-Keplerian angular momentum distribution along the disk. 
When more of the internal physics of the disk is included, 
this funtion may be determined. For some recent spherical 
accretion models r see Naraschi et. al. (1979) and Maraschi, 
Roasio and Treves (1981). 

An alternative possibility, for both the matter sur- 
rounding the hole and the method of energy release, is a disk 
that anchors a large-scale magnetic field, torquing and spin- 
ning down a rapidly rotating black hole (e.g. Blandford and 
Znajek 1977) . In this case the energy is supplied by the ro- 
tation of the hole. (At an earlier stage, this energy had to 
ultimately derive from energy made available in gravitational 
collapse.) Thermodynamically, this process can provide rela- 
tively low-entropy energy that is well suited for accelerat- 
ing electrons to relativistic velocities. 

C. Time Variability 

. . 3 

Time variability in the X-rays of ^yt 0 ~ S has 
now been observed in a number of QSOs and AGN. Some of these 
results are shown in the Table below, where the observed lu- 
minosity \-x is typically in the range 2-10 kev. 
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Table 1 


Source 

Lx 

AL^/L* 

At 


c At 


(erg/s) 


(s) 


(cm) 

Cen A 

Hi 

1x10 

0.25 

3 

7x10 


!H 

2x10 ! 

Hk 421 

1 mo” 

2.0 

1x10* 


3xl0 / * 

NGC 6814 

lxio” 

1.5 

2x10* 



OX 169 

1x10 

1.5 

6xl0 3 


, A \ 

1.8x10 < 

1 

3C 273 

1.7x10** 

0.1 

6xl0 3 


1.8x10’* j 

NGC 4151 

5x10 ¥Z j 

3.0 

lxlO 3 


3xl0’ 3 i 

1 


The observations for Cen A are from Delvaille et al. (1978) , 
for Mk 421 front Ricketts, Cooke and Pounds (1976) , for NGC 
6814 from Tennant et. al. (1981), for OX 169 and 3C 273 
from Tananbaum (1980) , and for NGC 4151 from Tananbaum et. 
al. (1978). a typical observed fluctuation is shown in Fig- 
ure 1 for Ox 169. 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 1: Einstein HRI observations of 
QSO OX 169 , from Tananbaum (1980) . 


The condition that the size of the emitting region dur- 
ing a fluctuation* H , be smaller than C A't can be written 



Since many of the fluctuations in Table 1 have Al^L * we 
can assume that a large i taction of the steady emission re- 
gion is involved in the observed fluctuations, . Note 
that, within the context of the black hole model, equation 
(10) and Table 1 are consistent for Mg -v, 0. 1 — I O . it 
is also important to note that if Mg***! • and the fluctua- 
tions are produced at a radius , the observed fluc- 
tuation timescales are shorter than an orbital period, cf. 
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equation (3) • Since gravitational energy cannot easily be 
released on a timescale £ , the energy in an outburst 
must have been stored. 

Some very general arguments may be given concerning the 
relation between L , AL , and At: . By assuming that the 
energy in an outburst is produced by material that is associ- 
ated with opacity, thereby increasing the light travel time 
across the emission region, Cavallo and Rees (1978) and Fabi- 
an and Rees (1979) have obtained the inequality 

AL ^ *. f c*e A-t/q t~‘ (£,) ■ 1 ") 

By requiring that the steady luminosity be less than the Ed- 
dington limit and using equation (10), Lightman, Giacconi, 
and Tananbaum (1978) have pointed out the inequality 

L < 10* a*' ^V'At . (.12.) 

Fortunately (for theorists), none of the observations in 
Table 1 violate the above inequalities, although 3C 273 and 
NGC 6814 push them. That these inequalities seem to have 
something to do with the actual data suggest both high effi- 
ciencies, £>0.l, and an emission region of dimensions deter- 
mined by the Schwarzschild radius of a black hole. 

Clearly, long-term observations with large area 
detectors and high time resolution are desireable for further 
testing the above inequalities and related considerations. 
The proposed Large Area Nodular Array of Reflectors (LAMAR) 
should be quite suitable for this task. 

For time resolution sufficient to study a signal of mag- 
nitude it may be possible to see large blobs of 
matter spiralling into the black hole. For a circular. 
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iwn-spi railing orbit, we expect to see angular frequencies 

\ I hj 

OJ- (6M /r 1 /,' in the Newtonian approximation, which can be 
written as . 

u'- u. ■ot*A,Y / * oVi 

where Co’o “ iM/ ? 0 is the angular frequency at radius f 0 and 
P is given in equation (3) . Now, for a radiating blob of 
matter that slowly spirals inward, the observed angular fre- 
quency will increase, seen as a low-Q quasi-periodicity, with 
decreasing period. To obtain an equation for this spiral, we 
may specify radius in terms of phase or azimuthal angle 0 , 
and then integrate the relation id - <1 tf/St. along with equa- 
tion (13a) to obtain 0 (t) . Signals would be modulated in 
proportion to cos (6. 

The type of spiral depends on the local, nongravitation- 
al physics, particularly the viscous stresses. As a simple 
example, the "Archimedes spiral" is of the form 

r- n, a- <f> \ 

2.7rn / 

where n is a free parameter, and yields 

0c-t) = ZKx\ f i- (i-fjeaS Y /S l Os i) 

Another spiral^w'Klch can be related to disk models in the li- 
terature^ derives from equation (5a), . If we 

assume that £ is a constant (as given by a number of the 
CL -disk models), then equation (5a) gives the logarithmic 
spiral 

r 0 e 

and yields 

0ft)' \ P'fi. 

3 r 3 a 

We mention that simple analytic solutions may also be obta- 




153 


ined when p is a general power law in radius ? rather than a 
constant. The above forms for ^ are only applicable several 
Schwarzschild radii away from the black hole; similar gener- 
al-relativistic expressions could be obtained closer to the 
hole. 

Further theoretical work needs to be done on the time 
variability to be expected by matter inhomogen.eities near the 
black hole. Accretion instabilites (Pringle? Rees? Pacholc- 
zyk 1973; Lightman and Eardley 1974; Shakura and Sunyaev 
1976) may contribute to the oberved variability. None of 
these possibilites has yet been modeled with sufficient deta- 
il? in the nonlinear regime? to make quantitative predic- 
tions. 

D. Gas Supply 

A clue to the environment of the power source might be 
gotten from study of the various sources of "fuel"? the re- 
quired conditions for these sources? and the different impli- 
cations for luminosity evolution. Some of the models that 
have been proposed for the source of gas are (a) stellar col- 
lisions in a dense star system (e.g. Spitzer and Saslaw 
1966) ? (b) tidal stripping of stars in a dense stellar system 
by a central massive black hole (Hills 1975) ? (c) infall of 
intergalactic gas (Gunn 1979)? (d) galaxy mergers (e.g. Roos 
1981) ? (e) intragalactic gas released by normal stellar evo- 
lutionary processes. One difficulty with models in which gas 
originates at large radii from the center is the necessity to 
dissipate a large amount of angular momentum. If the gas is 
provided by a dense stellar system? and if and 
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1>/J0 S'* , then it can be shown that the dynamics of the 
stellar system must be dominated by physical stellar colli- 
sions (Shields and Wheeler 1978; McMillan, Lightman, and 
Cohn 1981) . 

It seems clear that substantial luminosity evolution oc- 
curs in QSOs (e.g. Turner 1979); AGN may be the late time 
phase of the same evolution. In principle* the gas supply 
mechanism should have something to do with luminosity evolu- 
tion. The various mechanisms mentioned above, through their 
different dependences on the external environment of the QSO 
or AGN (in (c) and (d) ) , and on the internal parameters of 
the black hole and dense stellar system (in (a) and (b)), 
predict different functional forms for the evolving gas sup- 
ply rate WOt/ , and hence for the evolving luminc Ity L(t) . 
A comparison of theoretical models for the gas supply mechan- 
ism with a large data set of received fluxes and redshifts 
might be able to rule out or partly confirm some of the mo- 
dels and clarify the relevant physical conditions. For exam- 
ple, models (a) and (b) , with evolution of the stellar-system 
parameters included, yield late-time asymptotic laws of the 
form (McMillan, Lightman and Cohn 1981) 

M (-0 t' p , (.«•) 

where p®2 for mechanism (a) and p=l for mechanism (b) . 

The test of any theory of luminosity evolution is a sta- 
tistical problem, requiring a value for the cosmological par- 
ameter , some information about the distribution of ini- 
tial conditions or luminosities (the luminosity function) , 
and some information about the distribution of "turn-on 
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times" of the objects — in addition to an evolutionary law 
of the form of equation (14) for a single object. If addi- 
tional assumptions can be made about these distributions, 
they may not need to be completely specified a priori. For 
example. Turner (1979) assumes the luminosity function main- 
tains its shape in time, which, in our context, would allow 
determination of a value for p for every assumed value of 
. Unfortunately, Turner's assumption has little physical 
justification and seems to require some external clock that 
determines the physical conditions of each new quasar as it 
is born. In any case, a very dense X-ray Hubble diagram, with 
the fluxes and redshifts from many QSOs, will be necessary to 
solve the essentially statistical problem of luminosity 
evolution. X-rays seem particularly useful here, since QSOs 
may be readily identified by their X-ray emission, and this 
emission may constitute the majority of the total luminosity. 

h good determination of the local luminosity function 
would be quite useful in the above project. Figure 2 below 
shows a luminosity function calculated for a number of nearby 
objects. 
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o X-ray luminosity function of Seyfert jalaxies . 

• Local , optical luminosity function of QSO’s. 

• Space density of Seyfert nuclei f optical) . 

* Radio luminosity c u notion of E+SQ jalaxies (1.4 GHz). 

♦ Radio luminosity function of E+SQ cores (S GHz). 

Figure 2: Luminosity function of QSOs 
and AGN. The dashed line is fit from 
model (d) and additional assumptions. 

From Roos (1981) . 

The universality of the forms in Figure 2, especially above 
about L«3*10 , suggests some universal type of process that 

establishes the range of conditions for QSOs and AGN. 

III. RADIATION PROCESSES 
A. Types of Mechanisms 

It should be stated at the beginning that, while much of 
the emission in the radio is easily diagnosed as synchrotron 
radiation, there is little agreement on how the X-rays and 
gamma rays are produced. It is almost certain that the con- 
ditions in the X-ray emitting region are far from simple. A 
good fraction of the particles, especially those that ulti- 
mately form the observed jets, may be relativistic; there 
are accretion instabilities; shock waves and magnetic fields 
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abound; the emission may be part thermal and part nonther- 
mal. The radiation mechanisms we discuss can be divided into 
two categories: 1. Power-law electron models and 2. 

Non-power-law electron models. The second category can be 
further subdivided into thermal (random motions of electrons 
dominate emission) and nonthermal (bulk motions of eletrons 
dominate emission) processes. 

As will be shown below, the observed X-ray emission is 
well fit by a power law of the form 

I v « V s , 

where Xy is proportional to the energy per unit frequency V . 
1. Power-law electron models 

If the X-ray emission is synchrotron radiation (requir- 
ing electron Lorentz factors for equipartition mag- 

netic fields of 6 ) or relativistic inverse Compton, 

then a spectrum of the form of equation (15) can be achieved 
by a power law distribution of electrons of form 

dn./jy * y T tifo'i 

where 

For the observed values of the power-law index, S^0.7 ,see 
below, we require \r~2.H • A number of groups (Axford, 

Leer and Skadron 1977, Bell 1978, Blandford and Ostriker 
1978) have investigated first ^r^er Fermi mechanisms, whereby 
particles are accelerated into a power law distribution by 
the converging flows across a shock front. These mechanisms 
yield values f~2-3 , and have been applied to cosmic rays. 
It is also quite possible that such processes accelerate 
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electrons in the central regions of AGN and QSOs, where 
strong shocks may be present , and evidently produce an elec- 
tron index in the appropriate range for synchrotron or in- 
verse Compton radiation. If this is the nature of much of 
the radiation produced, then high X-ray polarization might be 
expected. Future X-ray plarization experiments may be quite 
important in clarifying this possibility. 

2. Non-power law electron models 

A power law spectrum can also be produced by the 
scattering of soft photons transmitted through a finite medi- 
um of hot electrons, even when these electron do not them- 
selves h*ve a power-law distribution. In this case, a power 
law results fcom the absence of an enerqy scale over a large 
range of energy. When the electrons have a thermal, nonrela- 
tivist ic distribution, then s may be related to the Thomson 
depth T.^ and temperature T (Shapiro, Lightman and Eardley 


1976, Katz 1976) L t+K 0 



£,= -! «•(.£. ity/* 

(I7«0 

where 

ys *f-T, 

(. I7L) 


T, 2 WT/md 2 . 

c ife') 


When the electrons have a thermal, relativistic distribution, 
then the corresponding result is (Pozdnyakov, Sobol and Suny- 


acv 1976) m h <l) 

s « r;; ci6 t, 1 ) u*) 

We mention that power laws may also be produced by the above 
process for non-thermal electrons, as long as their distribu- 
tion is sharply peaked at some "effective temperature.” Taka- 
hara (1980) , Takahara, Tsuruta and Ichimaru (1981) and Maras- 
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chi, Roasio and Treves (1981) have recently considered models 
using the above radiation mechanism. The soft photons for 
this process could be generated by cyclotron radiation (Eard- 
ley and Lightman 1976, Takahara, Tsuruta and Ichimaru 1981) , 
produced at perhaps /O* 0 Hz and absorbed up to 
about the 100th harmonic. Soft photons could also be pro-* 
duced by cold dense clouds in pressure equilibrium with the 
hot phase of the plasma. 

A problem with the Comptonization mechanisms above is 
that the spectral index depends on the parameters TT+t, and 
T , which might be expected to vary from one source to 
another, in contradiction with the narrow range of s actually 
observed (see D. below) . For a scattering model to work in 
this context, a "characteristic" geometry or structure is re- 
quired. For example, if the soft photons are externally in- 
cident on a region of hot, thermal electrons with a large 
scattering depth, then the reflected X-rays have a more univ- 
ersal spectrum. Lightman and Rybicki (1979a) have investi- 
gated this process and obtain the universal quasi-power law 
( for -V V? ) 

S - X ( tnHY 

2. Vo' J 

where V 0 is the typical frequency of the soft photon input. 
This power law extends up to photon energies and into 
the relativistic domain if the electrons have relativistic 
random motions. 

In the above scattering mechanisms, it is the random mo- 
tion of the electrons that dominates the emission. A closely 
related mechanism that depends on the bulk electron motion is 
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the scattering of intially soft photons by the converging 
flows across a shock front. This process has been investi- 
gated by Blandford and Payne (1981a, 1981b) and gives a high 
energy asymptotic result 

s — (n’-'/z; (n x +c) Citi 

where H is the Mach number of the shock (not to be confused 
with a mass) . This power law extends up to frequencies such 
that the photon momentum equals that of the electron. Thus, 
the bulk motions must be relativistic to produce photons of 
energies Vw > me** . Since characteristic bulk velocities 
are only of the order v*~»0.3c at r~tO , this scattering 
process may have difficulty in producing relativistic pho- 
tons. In any case. Equation (20) gives a value of s a little 
too sm.^11 and equation (21), foi a strong shock, M>",1, gives 
a value of s a little too large. It is possible that refine- 
ments of these models might be more satisfactory. 

Some general statements may be made. Any thermal pro- 
cess must result in a characteristic thermal (exponential) 
turnover in the spectrum, and possibly a Wien hump, at 
hV f* kT. A shock mechanism should have a cutoff in the 
spectrum at energies rnc.(Tv)^.^in almost all cases, no 
turnovers of any kind have yet been seen. It is important to 
extend X-rav observations up to higher energies, in search of 
such a turnover . In any case, it is clear that some turnover 
must occur, since the observed value of s^0.7 indicates an 
infinite total energy if this power law continues to arbi- 
trarily high frequencies. 

Changes of the spectrum during luminosity fluctuations 



may also provide some information about the emission mechan- 
ism. In any scattering model in which soft photon are scat- 
tered up to high energies r one expects the spectrum to harden 
during a fluctuation if the temperature remains constant. 
Rybicki and Lightman (1979b) and Payne (1980) have investi- 
gated the time-dependent versions of Comptonization by ther- 
mal electrons within this context. Some attempts have been 
made to model the observed (or absence of) spectral evolution 
in outbursts by such processes (e.g. Tennant et. al. 
1981) . If the Compton cooling time is much shorter than the 
energy input time r it may be the soft photon source rather 
than the temperature that remains constant during fluctua- 
tions. Guilbertr Ross and Fabian (1981) have investigated 
this possibility and find that the spectrum first hardens and 
then softens as the temperature decreases. Over a range of 
parameters with y£l • the time-averaged spectrum has a 
quasi power-law shape, with 0<s<l, and a thermal cutoff (with 
no Wien hump) at the initial electron temperature. It is 
possible that in many cases, steady state spectra are not re- 
levant and we should always be considering time averages of 
nonsteady processes. Better spectral resolution during 
luminosity fluctuations will be helpful here. 

There are various pocesses that may be involved in pro- 
ducing the electron distribution, if it is not thermal. In 
addition to the shock acceleration already mentioned, large 
scale magnetic and electric fields may be present (Blandford 
1976, Lovelace 1976, Blandford and Znajek 1977) . Cavaliere 
and Morrison (1980) have pointed out that the electrons must 



r . *3 

^ - f O’-' ' < V *- ! 

be reaccelerated in crossing the emission region. Defining a 
parameter q by 

q«/ probability \ ^/fractional energy loss 
Ipf scattering) v P er scattering 

= x. ^ i» J j jnv/jfrmcO • 

Here Df is the photon density 

Hy - L / (H'TTr^cUv) ? 

so that 

70 * (tj) &. 

In any source where q exceeds unity, reacceleration is neces- 
sary. Note that q is about equal to the Kepler time divided 
by the Compton cooling time. 

B. Thermalization Time Scales 

A principal uncertainty is whether any of the observed 
X-ray emission from QSOs and AGN is thermal. Such emission 
requires that the electrons be able to thermalize r via 
two-body collisions or perhaps collective plasma motions, on 
a time scale shorter than energy is deposited or radiated and 
shorter than pair creation when the latter is important. We 
will consider two-body thermalization below, keeping in mind 
that collective effects may always be faster. 

First we consider the ion thermalization time. The bulk 
of the released gravitational energy is probably intially de- 
posited into the ions. Taking a value of 25 for the Coulomb 
logarithm, the ion thermalization time, tjj , via ion-ion 
scattering is, cf. Spitzer (1962) 

t Vl = 5 *i0 ,7 s N*' ( 

Substituting equation (6) into equation (25a) , and taking T» 




io'*K 


). (20 


C Z2) 

C2-Z) 

(zO 
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10 K, we obtain 


t;; = Sato’s (j^) P £o.i) M 8 > 
ising equations (3) and (5b) , we obtain 

t,;/tr = I.C.XI0 2 (^)V (fl). 


(,2-Sl) 


U sO 


Here and elsewhere in this section, all quantities will be 
evaluated at f^lOv^ for simplicity, unless otherwise stated. 
Since energy is deposited into the ions on a timescale fc r , 
the last ratio must be less than unity for ion-ion collisions 
to thermalize the ions. Ions may also thermal ize by scatter- 
ing with the electrons. If we denote this timescale by 
•fc.V , then . x _*/- 

tii /tie v Cr>ip) ( N / ' ^ » 
where T-, and ~le denote the ion and electron temperature, res- 
pectively. (If the particles are not thermal , then tempera- 
tures should be replaced by mean energies.) For and 
Te. = K, equation (26) gives ^ 50 (n-tfO/W- 
Thus, a relatively low electron mean energy and the presence 
of many electron positron pairs may allow the ions to ther- 


malize. 


Next we consider electron 


thermal izat ion 


electron-electron scattering. For relativistic electrons, 

T. ^ I r this tilQ6SCdlGf t*c r IS 

t. e = Tj (n-HO' l 6:<r T jJ>.Ay l C 27 ) 

- *■**»** T ' ( ' (tb) p (f i ) n * • 


The electron cooling time, ir coo i * 
model-independent way: 


may be evaluated in a 
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t«i ^ |7tr3(o-tn4C2kT)/t. = T« (SiEt) (£) (/3e)\. (23 i 

We then obtain the ratio ( 

~ 0. 1 T* C~zf (£-J (3 ! (£[)*, «?) 

which must be less than unity for the electrons to thermal- 
ize. Equation (29) suggests that in some cases we might ex- 
pect thermal electrons and in some cases not. As long as the 
electrons are only marginally relativistic, the electrons 
should be able to thermal ize unless the luminosity is very 
much below the Eddington limit and there are few electron po- 
sitron pairs. It seems likely, however, that a fraction of 
the electrons will be accelerated up to very relativistic en- 
ergies by some of the processes mentioned above. In this 
case we would not expect the electrons to be able to sustain 
a thermal distribution. In any event, as is discussed below, 
the copious production of pairs prevents the temperature 

of an optically thin thermal medium from becoming very rela- 
tivistic. 

The time scale for electrons to receive energy from 
ions, tv, , satisfies 




Since electrons can cool relatively efficiently, faster than 
they receive energy from the ions, we may expect the electron 
mean energy to be much lower than that of the ions, as sug- 
gested by Shapiro, Lightman and Eardley (1976) . 


The conclusion is that there may well be two populations 
of electrons: a marginally relativistic, thermal population 

and a highly relativistic, nonthermal population. 
Unfortunately, it is difficult to calculate the relative pro- 
portions of these two populations at the present time. 
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C. Pair Effects 

With the high virial temperatures available, cf. 
equation (2) , much exceeding the electron rest mass, it seems 
likely that some sizeable fraction of the electrons will be 
at least marginally relativistic, and most detailed models in 
the literature (e.g. Shapiro, Lightman and Eardley 1976; 
Maraschi, P.oasio and Treves 1981; Takahara, Tsuruta and 
Ichimaru 1981) do indeed have this character. In such a si- 
tuation, the effects of electron-positron pair production 
must be incxided. When pairs are present, the distributions 
of radiation and pairs must be solved for self consistently, 
since photons produce pairs, pairs produce photons^ and all 
the cross sections are, in general, energy dependent. Some 
recent work in these processes has been done for steady plas- 
mas by Stoeger (1977) , Liang (1979) , Lightman and Band 
(1981) , Lightman (1981) , and Svensson (1981) , and in the con- 
text of a time-dependent cooling "fireball" by Cavallo and 
Rees (1978) . 

As an example, we consider production of pairs by the 
reaction 

er + X — » e x + +e" 9 

which may dominate when a large photon density is present and 
the electrons are energetic. This reaction gives a positron 
production rate 

CO^*r (ru*ru) 

(where off is the fine structure constant)^ only logarithmi- 
cally sensitive to the particle and photon distributions, as 
long as the threshold conaition 
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is strongly satisfied. The pair production time scale for 
this process, "fcpeY , is 

= r V/n + =• . 06 t r (3 C LTd 

where we have used equations (23) , (5b) and (3) , and the L 

refers to the luminosity in those photons above threshold. 
In sources that are not too far below the Eddington limit f we 
therefore expect this process to have ample time to produce 
pairs. Since pairs produce more pairs, there may be a non- 
steady behavior in which, for a short period of time (and 
with Ay varying rapidly perhaps), an overabundance of pairs 
is produced. 

It is also of interest to compare the time scale, 
-tr ttnn • for the pair annihilation process 

e — *■ zt 

to the cooling time. (In a steady state the creation time 
would equal the annihilation time.) If pairs cool before they 
annihilate, then we might expect to see a feature in the 
spectrum at the electron rest mass energy, as has been seen 
(redshifted by 20%) in the March 1979 gamma ray burst event. 
For QSOs at redshift z, this feature would be at the energy 
511 kev/(l+z). For electrons of Lorentz factor ^ , 

U>vV(27<r T cn.). <3z) 

In addi 'on to comparing equation (32) to equation (28) we 
can assume a specific cooling mechanism. For example, the 
time scale for synchrotron cooling, » is 

CKnc /W 3V t ) . (35) 

Then, using the expression for the gas and magnetic pres- 
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sures, p & , and pg , respectively, 

p G = L3?) 

Pb = BVb rc , 

we can obtain the ratio 

tvn/t artn = 32T~* ( P<&/Pg) • 

For equipartition magnetic fields, Pg ** Pg , pairs will ev- 
idently cool before they annihilate as long as / X 2 . The 
absence of an observed annihilation line may indicate the 
magnetic fields are below their equipartion value. When re- 
lativistic pairs annihilate before substantial cooling, the 
photons produced will at most yield a broad feature at 
h J Yfr. cP - , and even this may be dominated by Comptonized 
bremsstrahlung radiation (Lightman and Band 1981) . 

Recent work (Lightman 1981; Svensson 1981) has shown 
that there is a maximum temperature that can be achieved by a 
thermal, optically thin plasma in equilibrium, dependent on 
the value of the ion scattering depth, equation (7) . A maxi- 
mum temperature comes about because the pair annihilation 
rate decreases with increasing temperature and the pair pro- 
duction rate increases as the photon density increases (which 
increases with ~Lh ). This maximum (or, equivalently, the max- 
imum "0(4 possible for each value of T) is shown in Figure 3. 
Furthermore, there is generally an enormous range in luminos- 
ity over which the temperature remains in the narrow range 
O-l $7* £ I , with copious pair production at the high end and 
threshold effects at the low end serving as a thermostat. 
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Figure 3: Maximum ion scattering 
depth , Thomson depth, and photon 
to electron ratio that can be 
achieved at each temperature in 
an optically thin, steady thermal 
plasma. From Lightman (1981) . 

Earlier work by Bisnovatyi-Kogan, Zeldovich, and Sunyaev 
(1971), which did not include finite size effects and ne- 
glected photon processes, obtained the maximum temperature 
kT*41mc i . The magnitude of this maximum is set by the reci- 
procal of the fine structure constant. In the calculations 
resulting in Figure 3, only internal sources of photons were 
included (bremsstrahlung. double Compton, and annihilation). 
Any additional source of photons, e.g. a soft photon source, 
must lower the maximum T at each value of , since photons 
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produce but d v not destroy pairs in an optically thin plasma. 
Thus, the maxima shown in Figure 3 are upper limits. 
Observed spectra with power laws extending above 
lr\V > kirn** indicate that either the plasma is nonthermal 
or nonsteady. 

Another constraint on the mean energy and number density 
of the pairs arises from requiring that the total system of 
ions and pairs must be gravitationally confined. This leads 
to the inequality 

Y( n -~r-)^ t(%) ( tJ ■ w) 

When equation (37) is violated, magnetic confinement requires 
magnetic fields above equipartition strength. Equation (37) 
must be satisfied on average, but could be violated for short 
times or over small regions. 

We finally mention that photon-photon pair production 

may limit the presence of nard gamma rays. If is the 

energy of a soft photon and is the energy of a hard pho- 

ton, then the threshold condition for pair production is 

(.**.') / JuSlN >, 2 

v V- me 1 / ' 

The scattering depth to the above reaction is 

where is the number density of "field" photons in a 
bandwidth around Ve . Some constraints may be placed on 

sources for which strong gamma emission is seen, by requiring 
that the number density be sufficiently small to achieve 

, (Herterich 1974). For example, Fabian and Rees 

(1979) have shown that the observed gamma and X-ray emission 
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from 3C 273 (see Figure 4) require that the gamma rays be em- 

• 6 

itted in a region 1^710 cm , much larger than the X-ray emit- 
ting region. We mention, however, that the above constraint 
may be violated if the produced pairs themselves produce more 
gamma rays that are not substantially downgraded in energy. 
Since the various reactions that would redistribute the spec- 
trum, e.g. Compton scattering, are energy dependent, it is 
not clear how large can be before the gamma rays would 
certainly be cut off. This is an area in which more research 
is needed. 

D. Observed Spectra 

The composite spectrum of the QSO 3C 273 is shown in 


Figure 4 
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Figure 4: Spectrum of 3C 273. From Bradt (1980) 

The portion of this spectrum from 13 to 1L0 kev is well fit 
by a power law with s»0.67. This s^/ectrum is actually fairly 
characteristic of many AGN and QSOs. In Figure 5, we see a 
histogram of spectral ^ndicies that describe power laws seen 
in Seyfert galaxies. (The index oG in ?igure 5 is the same 
as the s given i;. iyuation (15).) 
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Figure 5: Histogram of spectral indicies observed 
for Seyfert galaxies, oc is the energy index. 

From Mushotzky (1981). 

The data, over an observed range 1-100 kev, are well fit with 
an s 0.69, and a srall scatter about this mean C~ ^0-/5" . 
Th*s result has been reported by Mushotzky (1981) and by Mac- 
cacaro, Perola and Elvis (1981) . The universality of these 
spectra is striking and suggests a universal type of emission 
mechanism, as discussed above. Emission above 1 Mev has also 
been seen in a u. other ob 3 ects, e.g. NGC 4151 (Schonfelder 



1978) . 




of AGN. For t' s purpose large area detectors like LAMAR 
should be useful. 

From all the a bove considerations, there is a clear need 
for observations at higher energies , going into the gamma 
rays. There may be a wonderland of interesting physics wait- 
ing to be learned at hv 'Z, mc* . The proposed satellite 
GRO should have such a capability. To summarize the needs 
for such observations: 

(1) With a spectral index of s»0.7, most of the total lu- 
minosity has yet to be observed! Look back at all the 
go-or-no-go formulae that contain an L. Until the observa- 
tions are pushed to sufficiently high energies to see a tur- 
nover in the spectrum, we only have lower limits on that L. 

(2) When a break, or turnover, in the spectrum is seen at 
high energies, we will have a much better idea of whether the 
emission is thermal or nonthermal. 

(3) As suggested by the discussion in C. above, a number 
of interesting effects associated with electron-positron 
pairs may be seen at energies Uv 2 • These features and 
the physical effects they represent may be completely hidden 
at lower energies. 

Observations of the spectra at low energies, especially 
near the soft X-ray absorption cutoff, 1 kev, are also 
important. Lawrence and Elvis (1981) have recently found 
good evidence for a correlation between total X-ray luminosi- 
ty and the ratio of hard ( hv- 2-10 kev) to soft 
( hv>0.5-4.5 kev) components. They interpret this result to 
indicate that the covering factor of the broad line emission 
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region over the continuum X-ray source is a mono toni cal ly de- 
creasing function of X-ray luminosity. Two possible explana- 
tions are that brighter sources, either through radiation 
pressure or hydrodynamic effects, are able to "blow away” the 
surrounding absorbing gas, or that the accumulation of ab- 
sorbing gas is associated with ageing (e.g. through stellar 
evolutionary processes) and older objects are in the declin- 
ing phase of their activity (see luminosity evolution in II 
D. above) . 

IV. BEANS AND JETS 

It is impossible to consider the subject of QSOs and AGN 
without mentioning the dramatic, well collimated streams of 
matter seen emanating from the centers of these objects. 
These "jets" have been detected in a variety of different ob- 
jects and in several different wavelengths. Recently, they 
have been seen in the X-ray band also (e.g. Schreier et. 
al. 1379, Teigelson et. al. 1981). The opening angle of a 
jet is typically Q £ 10° , the velocity of the matter 
300 V<C r and the jet lifetime 10 y*. See Rees 
(1980) and Begelman, Blandford and Rees (1981) for recent re- 
veiws of the physics of jets. 

Since jets are ubiquitous in QSOs and AGN, whatever pro- 
duces a jet should be a commonly occur ing phenomenon, not re- 
quiring special conditions. Furthermore, the jet "remembers” 
its directionality for a million years or more, so the funda- 
mental beaming mechanim must be associated with a longlived 
and stable axis of symmetry. Most researchers feel that net 
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angular momentum is a key ingredient in defining the required 
axis, but the detailed mechanism for accelerating matter 
along this axis is not yet understood. 

One of the earliest suggestions was the "twin-exhaust” 
hydrodynamic model of Blandford and Rees (1974) , in which a 
heat source at the center of a rotating gas cloud, in an 
external gravitating potential , creates a low density channel 
of gas that burrows a tunnel along the path of least resis- 
tance* the rotation axis. The first detailed bydrodynamical 
calculations to test the validity of this suggestion have re- 
cently been completed by Norman et. al. (1981), and some of 
their results are shown in Figure 6. 



Figure 6a: Density contours of 
the cloud jet structure at t» 
1*6. The solid line is the jet 
wall. The rotation axis is ver- 
tical. Units given in the text. 
From Norman et. al. (1981). 


Figure 6b: The velocity 
vector field for the jet 
in 6a. Note the central 
shock. From Norman et. 
al. (1981). 
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The unit of length here is GM/; , where H is the mass of the 
confining cloud and G. is the internal energy per unit mass 
of the hot, light gas introduced at the center. The unit of 
time is GM . The central cavity is at a radius R=0.15. 
The gas accelerates smoothly to supersonic velocity, termi- 
nating in a shock front at R*0.1, at which point collimation 
begins. These results indicate a cavity-nozzle structure 
more compact than that suggested by Blandford and Rees 
(1974). The calculations of Norman et. al. (1981) assume a 
flat bottomed gravitational potential, rather than a point 
source potential that would be generated by a central black 
hole. Future such calculations, perhaps more applicable to 
the AGN and QSO context, will yield one of the first reliable 
estimates for the amount of collimation that may be achieved 
by this mechanism. 

Another idea, suggested by Lynden-Bell (1978) , is that 
radiation pressure may accelerate matter along the interior 
funnel formed by a thick accretion disk (e.g. Jaroszynski, 
Abramowicz and Paczynski 1980) . Some recent kinematical stu- 
dies by Abramowicz and Piran (1980) suggest that the disk 
funnel must be extremely narrow for this method to yield ade- 
quate collimation, with a ratio of outer to inner disk radii 
of at least (O . Po calculations have yet been done to 
self consistently determine the parameters of the tunnel, or 
to test whether such a narrow funnel is stable in the pres- 
ence of radiation and matter. 

It is quite possible that none of the above mechanisms 
are able to achieve the observed collimation, and that much 



of the collimation is actually achieved far from the black 
hole, perhaps by large scale magnetic fields. Blandford and 
Payne (1981c) have recently investigated such a model, in 
which angular momentum is removed from an accretion disk by 
magnetic fields that extend to large radii. An outflow of 
matter along the rotation axis is then pcssiole. At large 
distances from the disk, the toroidal component of the mag- 
netic field may collimate the matter. 

V. CONCLUSIONS 

We conclude by recapping some of the important future 
observations that may clarify the physics of QSOs and AGN. 

1. Temporal Resolution, (a) The various theoretical 
constraints on L, L L, and At, when pushed to their limits, 
may provide our best handles on the efficiencies and 
emission-region sizes. (b) Spectral changes during fluctua- 
tions may help specify the radiation mechanism. (c) Signa- 
tures of spiralling matter or accretion instabilities may be 
revealed by high time resoltuion. 

2. Large Sample of QSO Spectra, Intensities, and Red- 
shifts. (a) It is important to learn whether a single univ- 
ersal spectrum describes QSOs as well as AGN. (b) With a 
large X-ray Hubble diagram for QSOs, we can begin to test lu- 
minosity evolution against theoretical models of the gas sup- 
ply mechanism. 

3. X-ray Polarization. Synchrotron radiation typically 
produces large polarization. Most other mechanisms produce a 
small polarization. 



4. High Energy Measurements V\V~ l-iO . (a) With 
an observed spectral index S^0.7 tor AGN, we now have only a 
lower bound on the total luminosity of most of these objects, 
(b) The shape of the turnover in the spectrum, at high ener- 
gies, should tell us about the emission mechanism. (c) A 
good deal of physics associated with pair production (maximum 
temperatures of thermal, steady plasmas; annihilation fea- 
tures) can only be observed at energies l\V>.SMcV. 

To end with a note of caution: Most of the ideas have 
been presented as if there were one universal framework for 
all QSOs and AGN, one basic scenario of gas accretion onto a 
massive black hole, one dominant radiation process, one me- 
chanism for producing and collimating the jets. I think 
sometimes, especially in astronomy, we tend to use Occam's 
razor too much. We should be receptive tc the possibility 
that a number of different processes and mechanisms may be 
operating in different objects, or perhaps going on simul- 
taneously in the same object — even though the net results 
of luminosities, spectra, and variabilities can be squeezed 
into a common framework. 
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